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ABSTRACT 



The nature of the asymmetry that gives rise to Type I X-ray burst oscillations on 
, accreting neutron stars remains a matter of debate. Of particular interest is whether 

^ ' the burst oscillation mechanism differs between the bursting millisecond pulsars and 

the non-pulsing systems. One means to diagnose this is to study the energy dependence 
of the burst oscillations: here we present an analysis of oscillations from 28 bursts 
observed during the 2003 outburst of the accreting millisecond pulsar XTE J1814- 
338. We find that the fractional amplitude of the burst oscillations falls with energy, 
\Q , in contrast to the behaviour found by Muno et al. (2003) in the burst oscillations 

' from a set of non-pulsing systems. The drop with energy mirrors that seen in the 

^ c*| accretion-powered pulsations; in this respect XTE J1814-338 behaves like the other 

accreting millisecond pulsars. The burst oscillations show no evidence for either hard 
or soft lags, in contrast to the persistent pulsations, which show soft lags of up to 50 
!_i ' (is. The fall in amplitude with energy is inconsistent with current surface mode and 

£^ \ simple hot spot models of burst oscillations. We discuss improvements to the models 

and uncertainties in the physics that might resolve these issues. 



Key words: binaries: general — stars:individual (XTE J1814-338) — stars:neutron- 
■ stars:rotation — X-rays:binaries — X-rays:bursts 



1 INTRODUCTION 

The neutron star Low Mass X-ray Binaries (LMXBs) are neutron stars accreting from low-mass companions that over-fill their 
Roche lobes. Of these systems some, the accreting millisecond pulsars, are also observed as X-ray pulsars. Of the « 160 known 
non-pulsing LMXBs, nearly half exhibit Type I X-ray bursts due to unstable nuclear burning of accreting material. High 
frequency modulations of the X-ray lightc urve, known as burst osc i llation s, have now been detected in the bursts from twelve 
of these systems (see the recent review bv lStrohmaver fc Bildsten IfeoOft h Of the seven known accreting millisecond pulsars, 
three exhibit bursts. Two of these systems, SAX J1808.4-3658 (hereafter J1808) and XTE J1814-338 (hereafter J1814), have 
burst oscillations. 

For the pulsars, the burst oscillation frequency is at or very close to the known spin frequency of the star Jchakrabartv et alJ 
2003: Istrohmaver et al. 1120031) . This suggests that the rotation modulates an asymmetry on the burning surface that is near- 
stationary in the corotating frame. For the non-pulsars the situation is less clear, because there is no independent measure 
of the spin frequency. The fact th at the frequency is highly stable from burst to burst, however ((Strohmaver et, al"1ll998l : 
iMuno. Ozel fc Chakrabartv l20o3h implies that there is at least a strong dependence on stellar spin. The detection of highly 
coherent oscillatio ns lasting several hundred seconds during a superburst from the LMXB 4U 1636-356 adds further support 
to this hypothesis l)Strohmaver fc Markwardt Il2002afl . 

The precise nature of the brightness asymmetry is as yet unclear. The simplest possibility is uneven distribution of fuel, 
leading to patchy burning. For the pulsars we know that initial fuel distribution is asymmetric, but the fuel could spread rapidly 
over the surface. For the non-pulsars there is no evidence of asymmetric fue l deposition, although there are several mechanisms 
that might render this process (and the associated pulsations) undetectab le iBjainer^i_fcJ_ ! anA ll987tlKvlafi~ fc Klimmis Il987t 
IWood. Ftaclas fc Kearney Tl98 8; Mcszaros. Riffcrt fc Berthiaume ll988llCummmg. Zweibel fc Bildsten l200ll : lTitarchuk. Cui fc Wood I 
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2002). Alternate mechanisms that do not rely on asymm etric fuel distribution include the development of vortices driven by 
the Coriolis forc e iSpitkovskv, Levin fc Ushomirskyj2002 |) . or global modes in the burning su rface layers iMcDermott fc Taaml 
Il987t ILeelbooi IrlevllbooilLee fc StrohmaveJbooa ICumming 1120051: IPiro fc Bildstenll2005l) . 

One means of diagnosing the nature of the asymmetry is to study the energy dependence of the burst oscillations. One 
can study both the dependence of amplitude on energy, and the relative phases of the waveforms as a function of energy. 
A "so ft lag" implies that the soft (low energy) pulse arrives later in phase than the hard (high energy) pulse, and vice 
versa. iMuno. Ozel fc Chakrabartv I <l2003l) examined the energy dependence expected for a simple hot spot model, generating 
simulated Rossi X-ray Timing Explorer (RXTE) Proportional Counter Array (PCA) lightcurves for neutron stars with a 
circular uniform temperature hot spot on a cooler background. The temperature contrast between the spot and the rest of 
the star, coupled with Doppler effects due to the star's rotation, leads naturally to a rise in fractional amplitude with energy 
in the 3-20 keV band. The hot spot model also predicts soft la gs. Models where the burst oscillat ions are caused by surface 
modes also predict a rise in fract ional amplitude with ene rgy lHevlll2005l : I Piro fc Bildsten"ll2006l) . In most cases the mode 
models predict soft lags, although iLee fc Strohmaverl (120051) found hard lags for a restricted subset of possible geometries. 

The most comprehensive study to date of burst oscillations from the no n-pulsing LMXBs, using RXTE PCA data, found 
that fractional amplitude rises with energy over the 3-20 keV energy band (|Muno. Ozel fc Chakrabartv II2OO3} ) . The authors 
show that the observations are compatible with the simple hot spot model if the bulk of the stellar surface (away from 
the hot spot) emits in the lower part of the PCA energy band. The data are also compatible with surface mode models 
iPiro fc Bildsten l Eool). provided that the required amplitudes can be excited. 

Muno et al. also examined whether there were phase lags between soft and hard photons. For most bursts the data were 
consistent with there being no phase lag; for 13 of 54 bursts studied the data were inconsistent at the 90% confidence level 
with being constant. In these bursts the trend was for the hard photons to lag the soft photons by up to 0.12 cycles, whereas 
both simple hot spot and mode models predict soft lags. The authors considered several mechanisms that might reverse a 
soft lag. They concluded t hat inverse C ompton scattering of soft photons to higher energies by a hot corona of electrons was 
the most likely candidate iMiller Ill995l) . The observational implications of such a corona, however, and its behaviour during 
a burst, remain to be worked out in detail. 

To date there has been no study reporting the energy dependence of the burst oscillations of the accreting millisecond 
pulsars. Given the possibility that the burst oscillation mechanism may differ from that operating in the non-pulsing sources, 
this is clearly an omission. In this paper we attempt to rectify this by presenting an analysis of the energy dependence of the 
burst oscillations of J1814. 

The source was first detected in outburst on 2003 June 3 by the RXTE Galactic bulge monitoring c ampaign. Its status 
as a millisecond pulsar was confirmed by a longer observation on June 5 iMarkwardt fc Swank 1 120031) . The pulsar has a 
spin freque ncy of 314.36 Hz, resides in a binary w ith orbital period 4.275 hr, and has a minimum companion mass of 
~ O.15M0 ijlVIarkwardt. Strohmaver fc Swank Il2003f) . The binary orbit is the widest of the seven known accreting millisecond 
pulsars. J1814 remained in outburst until mid- July 200 3, and in this time 28 Type I X-ray bursts were observed, all with 
detectable burst oscillations at the spin frequency iStrohmaver et al~ll2003( l. The bursts show significant harmonic content, 
wh ich could be used to constrain the equatio n of state and the system geometry fehattacharvva et al. 120051) . A detailed study 
bv lWatts. Strohmaver fc Markwardt I pOORl ) found that there is no evidence for fractional amplitude variation or frequency 
shifts in any of the bursts apart from the one burst that appears to show photospheric radius expansion. For most of the 
bursts, fractional amplitude is consistent with that of the persistent pulsations, although there is a small population of bursts 
with amplitudes that are substantially lower. For the first harmonic, substantial differences between the burst and accretion- 
powered oscillations indicate that hot spot geometry is not the only factor giving rise to harmo nic content in the latter . 
There are however no detectable phase shifts between the burst and accretion-powered pulsations l|Strohmaver et al~ll2003h . 
suggesting a model where the presence of the magnetic field somehow leads to a temperature asymmetry centered on the 
polar cap when the burst ignites (by allowing additional fuel build-up, for example). In this complementary study we test 
the models further by examining the energy dependence of the burst oscillations, comparing the behaviour to that of the 
accretion-powered pulsations. Section|5|details our method of analysis and our results. In Section[3]we discuss our findings in 
the light of current burst oscillation models, and outline areas for future theoretical study. 



2 DATA ANALYSIS 
2.1 Methodology 

Almost all of the PCA data for the outburst is event mode data with 125 fis time resolution and 64 energy bins covering 
the range 2 - 120 keV. The exception is the first burst, which was recorded in GoodXenon mode, which has higher time and 
energy resolution (lfis and 256 energy bins). Event mode data overruns, which are often seen in the bursts of brighter sources, 
were not seen in any of the J1814 bursts. The da ta were barycentered prior to analysis, using the JPL DE405 ephemeris and 
the source position determined from PCA scans IIMarkwardt fc S wank 1 120031: iKrauss et al"1l2005l) . 
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Computation of fractional am plitudes is done in two ways. We start by using the Z^ statistic feuccheri et al. Ill983l : 
Istrohmaver fc Markwardt 112002'bT) . This measure is very similar to the standard power spectrum computed from a Fourier 
transform, but does not require that the event data be binned. It is defined as 



^cosfc^j +|^ smfc< fej (!) 

where n is the number of harmonics summed (we use n — 1 throughout this paper), TV is the total number of photons, and j 
is an index applied to each photon. The phase 4>j calculated for each photon is 
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where u(t) is the frequency model and tj is the arrival time of the photon relative to some reference time. We use a frequency 
model in which the intrinsic spin rate is modified by orbital Doppler shifts using the best fit orbital ephemeris. Given a 
measured Z%, which we call Z m , the probability of the true signal power lying between and Z s is given by: 



f(Z s : Z m ) = exp 
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jGroth Ill97 j IVaughan et al~lll994h . We take the best estimate for Z s to be that for which f(Z s : Z m ) = 0.5. Given this Z s 
the RMS fractional amplitude r is then given by 



1/2 



N 



N-N b J (4) 
where Nb is the number of photons due to background accumulated in the energy band of interest during the observation 
period. Nf, is estimated using the standard FTOOLS routine pcabackest and the PCA background models. 

The second method used is to generate a folded pulse profile (using the frequency model) and then fit a sinusoidal 
model with as many harmonics as necessary. We use this method to check the fractional amplitudes computed using the Z„ 
statistic. Pulse profile fitting also allows us to calculate phase shifts between the different energy bands, something that is 
not possible using the Zt statistic alone . A more extensive discussion of both of these methods is given in Section 2.1 of 
IWatts. Strohmaver k. Markwardtl fe005h . 

The other issue to consider when computing fractional amplitudes for pulsar burst oscillations is that the accretion process 
may continue during the bursts. If this is the case then the measured fractional amplitude will contain contributions from 
both the burst process and the accretion process: 



r — 



f bur 



N. 



(5) 



iVbur and iVacc are the number of source photons arising from the burst and accretion processes respectively, with r^ UI and 
?~acc being the fractional amplitudes of the two different processes. The total number of source photons is N a = N^ur + A^cc- 
Since we are not always in the regime where Abur S> A^cc we will need to estimate 
rbur differs substantially from that measured. We do this in Section f2. 21 



and r acc in order to check whether 



2.2 Accretion-powered pulsations 

We start by folding together data from the whole outburst to generate a high resolution plot of amplitude and phase lags 
against energy for the accretion-powered pulsations. Figure Q shows the behaviour of amplitude against energy, for both the 
fundamental and first harmonic. There is a clear drop in amplitude with energy for the fundamental, of about 2% RMS over 
the 2-20 keV band. The amplitude of the first harmonic, by contrast, varies little with energy. Figure [5] shows the phase shifts 
between the different energy bands for the fundamental. Soft lags of up to 50 [is (« 0.015 cycles) develop over the range 2-7 
keV, with the lags levelling off at higher energies. The same behaviour is seen in the phase shifts computed using the first 
harmonic. 

This behaviour mirrors that seen in the accretion-powe red pulsations of the ot her millisecond pul sars: ,11808 |C ui 1 Mor gan fc Titarchuk 
ll99Sl: lOierhriski. Done fc Barret Il2002h . XTE J0929-314 JOallowav et al. Il2002h. XTE J1 807-294 jKirsch et al. 112004^ . IOR 
,100291+5934 jGallowav et al. Iboosir and XTE ,I1751-305 l|Gierlihski fc PoutanenlboOBh all show a drop of fractional am- 
plitude with energy within the PCA energy band. The other sources also show soft lags, although at only 50 /is the J181 4 
lags are smaller than those measured for the other systems. Detailed spectral modelling for J1808 bv lGilfanov et al. Nl998h . 
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Figure 1. Dependence of fractional amplitude on energy for the accretion-powered pulsations, computed from a folded profile using 
data from the whole outburst. The fundamental (top panel) shows a clear drop in amplitude with energy. The amplitude of the first 
harmonic, by contrast (lower panel) shows little change with energy. 



lOierliriski. Done fc Barret I l)2002h and IPoutanen fc Gierliriski I ( 20031) suggests that the amplitude drop and lags can be ex- 
plained by the presence of a hard Comptonized component due to boundary layer emission together with a softer blackbody 
component due to a hot spo t at the footpoint of the accretion chann el. Similar detailed modelling for J1814, which has a 
harder spectrum than J1808 jMarkwardt. Strohmayer fc Swank II2OO3I) . has yet to be done. 

In order to correct for the accretion contribution recorded during the bursts, we also need to know whether the properties of 
the persistent pulsations vary over the course of the outburst. Figure |3] shows the change in the countrate (excluding bursts) 
for the three energy bands 2-5 keV, 5-10 keV, and 10-20 keV. The evolution of fractional amplitude (at the fundamental 
frequency) is shown in Figure 2] The trend of amplitude dropping with energy is apparent over the whole outburst, although 
the magnitude of the drop does vary between 1% and 3% RMS over the outburst. The phase shifts between the different energy 
bands are shown in Figure Again the soft lags persist throughout, but get noticeably larger at the end of the outburst, 
when the accretion rate drops. It would be interesting to see if such a variation in phase lag with accretion rate is detected 
in any of the other millisecond pulsars. 



2.3 Burst oscillations 

The proportion of photons in the different energy bands for each of the bursts is shown in Figure [(J For most of the bursts 
the proportions are very similar. For the six faintest bursts, the proportions in the high er energy bands drop, as might be 
expected given that these bursts have lower peak temperatures l|Strohmaver et al. ^ 2003 ), Similarly, for the brightest burst, 
the proportion in the highest energy band rises. We define the duration of the burst as to be the time for which the countrate 
in a given energy band exceeds the ambient countrate by a factor of 1.5. Using this definition, burst duration is shorter at 
higher energies. 

We start by considering each of the bursts separately, computing fractional amplitudes and phase shifts for the three 
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Figure 2. Phase shifts compared to the phase for the lowest energy band (2-3.3 keV) for the accretion-powered pulsations. Negative 
values indicate that the pulse profile in the energy band being studied leads the profile in the lowest energy band, with 0.01 cycles 
corresponding to 31.8/is. The phase shifts were computed by folding together data from the whole outburst. 



energy bands 2-5 keV, 5-10 keV and 10-20 keV. The fractional amplitudes (at the fundamental frequency) for all of the 
bursts are shown in Figured and s umm arized in columns 2-4 of Table Q (for more general data on each burst see Table 1 of 
IWatts. Strohmaver fc MarkwardTI i2005t0 . More detailed plots for four example bursts are shown in Figure [H] We find that 
burst amplitude is lower in the 10-20 keV band than in the 2-5 keV band for 25 of the 28 bursts. Treating each burst as an 
independent test, we find that we can rule out the hypothesis that fractional amplitude is constant with energy at a level 
greater than 3a. The hypothesis that fractional amplitude rises with energy can be ruled out with an even higher degree of 
confidence. 

The phase shifts between the different energy bands for each burst are shown in Figure H There is no evidence for either 
hard or soft lags; the data are compatible with there being no phase shift. 

In order to verify these results, we folded together data from all of the bursts to generate a combined profile. We find a 
drop of fractional amplitude with energy of « 3% RMS over the 2-20 keV band, a steeper drop than that seen in the persistent 
pulsations (Figure [TTJ1 . The phase shifts of the combined profile, for two energy binnings, are shown in Figure [TTI The data 
are consistent with there being no phase shift, and we can rule out at the 3cr level hard or soft lags of greater than 0.015 cycles 
over the 2-20 keV band. We would have been sensitive to lags of the magnitude seen in the accretion-powered pulsations were 
such lags present. 



3 DISCUSSION 

The burst oscillations of J1814 exhibit a drop in fractional amplitude with energy of a few percent RMS across the 2-20 keV 
band. In addition the bursts show no evidence for phase shifts between soft and hard photons; we can rule out at the 3er level 
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Figure 3. Daily average countrate (excluding bursts) in different energy bands, corrected for background. Stars: 2-5 keV. Diamonds: 
5-10 keV. Triangles: 10-20 keV. 



shifts greater than 0.015 cycles over the 2-20 keV band. In this section we discuss these results in the context of the various 
burst oscillation models. 



3.1 Consistency with existing burst oscillation models 

We st a rt by considering models whe re surface modes cause the brightness asymmetry. The mode models discussed bv iHevll 
j2005h : Eeefc StrohmaveJ J2005h and lPiro fc Bildsten I &006) all predict a rise of amplitude with energy, although the gradient 
of the rise depends on various neutron star parameters. Since we can rule out a rise, or even a constant amplitude, with a high 
degree of confidence, the amplitude data are i nconsistent with these mo dels. The phase lag data do not provide any additional 
constraint, since the models of (for example') iLee fc Strohmaverl 1)20051 1 predict anything from hard lags of 0.04 cycles to soft 
lags of 0.1 cycle depending on the neutron star properties. Although we can rule out lags greater than 0.015 cycles there are 
still models that would fit the phase lag data. 

The other class of models that have been studied in any detail are hot spot models. These models are purely phenomeno- 
logical, in that the physics that might give rise to a hot spot remains as yet unclear. Like the mode models, the simple 
one-temperature hot spot models studied bv lMuno. Ozel fc Chakrabartv I <l2003l) also predict a rise of amplitude with energy, 
although again the magnitude of the rise depends on the neutron star properties. The models studied by Muno et al. also 
predict soft lags of least 0.02 cycles over the RXTE energy band. Any phase lags present in the burst oscillations of J1814 are 
smaller than this limit. 

We are therefore left with several options. One possibility is that the burst oscillation properties are consistent with 
existing mode/hot spot models but are masked by the presence of accretion-powered pulsations. We examine this in more 
detail in Section [3.21 and show that the accretion process would have to be disrupted substantially by the bursts to explain 
the observations. A second option is that existing models are not adequate to explain the observations. In Section 13.31 we 
explore whether a more complex surface temperature gradient might lead to consistency with the observations. 
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Figure 4. Fractional amplitudes for the accretion-powered pulsations for different energy bands, showing the variations over the outburst. 
Stars: 2-5 keV. Diamonds: 5-10 keV. Triangles: 10-20 keV. Each point is computed by folding data for a 5 day period starting on the 
day indicated. The drop in fractional amplitude with energy persists throughout the outburst, although the magnitude of the change 
does vary, particularly at the point where the accretion rate starts to drop. 

3.2 Masking by the accretion pulsations 

Unless the accretion process is completely disrupted by the burst, the measured fractional amplitude will contain contributions 
from b oth accretion and burst p ulsations. Separating the two contributions is not trivial because in this source they are phase- 
locked Jstrohmaver et al. 120031) . The persistent pulsations of J1814 do indeed show a drop in fractional amplitude with energy 
of 1-3% RMS over the 2-20 keV energy band (during the main portion of the outburst), so masking is not inconceivable. To 
proceed we need to estimate the magnitude of the accretion contribution to Equation 

We start by making the simple assumption that both the accretion rate and the amplitude/energy relation of the accretion- 
powered pulsations remain in their pre-burst state, as presented in Section f2. 21 during a burst. Columns 5-7 of Table^give an 
estimate of the ratio of accretion to burst photons over the duration of the burst, computed from the daily average accretion 
rate at the time of the burst. In Columns 8-10 we combine this estimate with the daily average fractional amplitude for the 
accretion-powered pulsations to estimate rb ur , using Equation 0. Although the accretion correction does make some small 
quantitative changes, it does not change the qualitative behaviour with energy for any of the bursts. Unless the accretion 
contribution differs dramatically at the time of the bursts from that observed normally, the burst fractional amplitude still 
falls with energy. So let us now consider the magnitude of changes to the accretion pulsations, compared to the properties 
measured immediately prior to each bursts, that would be necessary to mask burst pulsations whose fractional amplitude 
rises with energy. 

If the accretion rate remains at the pre-burst level, the fractional amplitude of the accretion pulsations must fall more 
steeply with energy during the burst than it does normally. Either the amplitude at low energies must rise, or the amplitude 
at higher energies fall, or both. For most of the bursts the ratio of accretion to burst photons is less than 0.3 (Table 0. The 
fractional amplitude would have to change by several % RMS to have the desired effect. Whether the burst process could lead 



8 Watts & Strohmayer 



<u -0.01 




d -0.03 

CL 



at aAaA a 



A 

A A a 



' 



20 30 
Time (Days) - MJD 52790 



40 



30 



Figure 5. Phase shifts compared to the phase for the 2-5 keV energy band for the accretion-powered pulsations, showing the evolution 
over the outburst. Top: 5-10 keV. Bottom: 10-20 keV. Negative values indicate that the pulse profile in the energy band being studied 
leads the profile in the 2-5 keV band, with 0.01 cycles corresponding to 31.8/is. Each point is computed by folding data for a 5 day period 
starting on the day indicated. 



to this level of disruption is not known, but it seems unlikely. A smaller change in fractional amplitude would be required 
if the accretion rate rose during a burst, boosting the ratio of accretion to burst photons. The effect of bursts on accretion 
rate is however unclear. Naively one might expect the radiation pr essure to hinder the accretion process. But the radiation 
may also remove angular momentum, increasing the accretion rate iWalker fc Meszaros 1 ll98St IWalker Ill99i iMiller fc Lamb I 
1993; IMiller lfl99i4 iBallantvne fc Everett 1 120051 . A substantial rise in accretion rate seems however unlikely, since the burst 
spectra are well-fit by a single temperature blackbody, and there is no spectral evidence for an enhanced accretion component. 

The one-temperature hot spot models, and many of the mode models, predict burst oscillations with detectable soft lags. 
This could only be masked by the accretion process if the accretion-powered pulsations exhibit hard lags during bursts. This 
would be a reversal of the normal behaviour, requiring substantial disruption of the accretion process by the burst. Whether 
or not this is plausible is a topic for future study. 



3.3 Hot spot with temperature gradient 

The evidence presented in Istrohmaver et al~l J2003) and IWatts. Strohmaver fc MarkwardTI i2005h suggests that the most 
promising model for the burst oscillations of J1814 is a hot spot that develops at the magnetic footpoint, perhaps due to a 
disparity in fuel build-up in this region. We will therefore foc us on the hot spot model in more d etail, and ask whether some 
simple alterations to the one-temperature hot spot models of iMuno. Ozel fc Chakrabartv I (I2003F) might enable us to explain 
the energy-dependences reported in this paper. 

Although most modelling to date has focused on uniform, single temperature hot spots, several processes involved in burst 
generation would seem to favour some degree of temperature variation across a hot spot. For example, we expect that burning 
must spread in some fashion such that regions which ignite first will have a longer time to cool than portions of a hot spot 



XTE J 1814-338 burst oscillation energy dependence 9 



0.6 



0.4 



o 



u 

o 



c 

o 

o 

JZ 



0.2 



0.0 



o 



X 



o 

o o 



X 



o 



o 

o o 



X 



o 

X 



* * X a X * X X ... * * * 



X 



X 



a aA a a a a a a a aA a a a a a a Aa 



A 



A 



A A A 



A 



A 



10 15 

Burst number 



20 



25 



30 



Figure 6. Proportion of photons in the different energy bands during each burst (corrected for background and an estimate of the 
accretion-related emission). Proportions are roughly constant, although for the faintest bursts (1, 6, 8, 9, 14 and 19, see Figure 4 of Watts 
et al. 2005) the proportion of photons in the highest energy bands drops. For the final, brightest burst, the proportion in the highest 
energy band is higher. Stars: 2-5 keV. Diamonds: 5-10 keV. Triangles: 10-20 kcV. 



which ignite later. The time asymmetry introduced by spreading can thus lead to a temperature asymmetry. Alternatively, if 
accretion is substantially mediated by the magnetic field, then the quantity of fuel could vary in a systematic fashion in the 
vicinity of the magnetic footpoint. Such a fuel asymmetry could also impart a temperature gradient within the hot spot. 

In order to explore the impact of a temperature gradient we have carried out model calculations using circular hot spots 
with a temperature profile across the spot. It is not our intention to fit detailed models to the data, only to examine whether 
a temperature gradient might allow us to fit the ampli t ude-e nergy relation. We use the model for emission of photons from a 
rotating neutron star surface described bv lStrohmaver I (120041) . The model includes bending of photon paths in a Schwarzschild 
geometry as well as relativistic beaming and gravitational red-shifts. Each surface element within the hot spot is assumed to 
emit, locally, a black body spectrum. The temperature distribution is azimuthally symmetric around the center of the hot 
spot, and is assumed to vary with the angle, a, measured from the central axis of the spot. For the purposes of this exploratory 
calculation we use a simple linear variation of the temperature with angle, kT = fcTo(l + Aa/ao), where kTo, A, a and ao 
are the central temperature, fractional temperature change, half-angle of the emitting point (measured from the spot centre) 
and maximum half-angular size of the hot spot, respectively. We compute a model photon energy spectrum at infinity, and 
then apply an RXTE/PCA response matrix to obtain count rates versus PC A energy channel. We then compute amplitudes 
as a function of energy channel that can be compared with the real data. 

Two example models are shown in Figure IT21 We have selected a geometry that gives an amplitude close to the observed 
values. As expected, and as shown by other authors, a uniform spot emitting a typical burst spectrum produces an increasing 
amplitude with energy in the PCA bandpass. Since in general a smaller hot spot produces a larger modulation amplitude, 
one would expect that in order to produce a drop in amplitude with increasing energy the temperature should increase away 
from the center of the spot. That is, the edge should be hotter than the center. This expectation is borne out in our modeling. 
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Table 1. Summary of burst properties. Photon ratios (columns 5-7) and accretion-corrected amplitudes (columns 8-10) 
are calculated assuming accretion rate and pulsation properties remain at pre-burst levels. No error bars are given for 
columns 8-10 because of the uncertainties inherent in this assumption: the error bars in columns 2-4 should be taken as a 

lower limit. 
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In the example shown, adding a fractional temperature change of ~ 10% produces a drop in amplitude with energy of similar 
magnitude to that observed. 

We do not claim that the above geometry is unique in fit ting this aspect of the data, a nd it is not clear that it would survive 
other constraints such as pulse profile fitting (see for example lBhattacharvva et al. I J200^)). However, it does demonstrate that 
temperature gradients merit further study. How such systematic variations in temperature might arise is not clear, although 
we note that an asymmetry associated with spreading from the center of a hot spot would tend to produce a hotter edge (as 
the center has longer to cool). Constraints from theory on the magnitude and distribution of likely temperature variations 
would be extremely helpful. 



4 CONCLUSIONS 

In this paper we have presented an analysis of the energy dependence of burst oscillations from the accreting millisecond 
pulsar J1814. Our results are intriguing. The fractional amplitude of the pulsations falls with energy across the 2-20 keV 
RXTE band; we are able to rule out the amplitude being constant, or rising with energy, at a level greater than 3a. This 
mirrors the behaviour seen in the persistent accretion-powered pu lsations, but differs from that seen for the burst oscillations 
of the non-pulsing LMXBs bv lfVIuno. Ozel fc Chakrabartv I l)2003h . The persistent pulsations, like those of the other accreting 
millisecond pulsars, show soft lags of 0.015 cycles over the 2-20 keV band; for the burst oscillations we can rule out lags or 
leads of this magnitude at the 3a level. 
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Figure 7. Burst RMS fractional amplitude (%) for different energy bands, at the fundamental frequency of the pulsar. Top: 2-5 keV. 
Middle: 5-10 keV. Bottom: 10-20 keV. The scales on each plot are the same to make clear the general drop in fractional amplitude with 
energy. 



The energy dependence of J1814's burst oscillations diffe rs from that reporte d for LMXB burst oscillations in the most 
comprehensive study of these objects to date jMuno. Ozel fc Chakrabartv II2OO3T) . It joins a list of properties that seem to 
differ between the two classes of objects. For the LMXBs, burst oscillations ar e only seen at high accretion rates; for th e pulsars 



for example iMuno et al. 


J200ll)'): for the pulsars frequency vari- 


dChakrabartv et alJ2003l: 


Watts. Strohmaver & Markwardt 12005): 



Bhattach arvva fc Strohm avcr 2006). The LMXB burst oscillations also show amplitude variations ( M uno l Ozel fc Chakrabartv 
2Tffl2rh for J 1814 no such variations are seen except during PRE llWatts. Strohmaver fc Markwardt! 200fir Finallv. the pulsar 
burst o scillations have detectable harmonic llChakrabartv et al.l2003l:IStrohmaver et al t 2003 ) , whereas the LMXB oscillations 
do not iMuno. Ozel fc Chakrabartv I J2002h . although see lBhattacharvva"etHlrt^T2o7)5^ . which presents evidence for harmonic 
content in the rising phase of bursts from 4U 1636-536). 

The differing properties suggest that the mechanism responsible for the burst oscillations may differ between the two 
classes of objects. Further study, however, is clearly required. There is now a much larger sample of LMXB burst oscillations in 
the RXTE archive, and it would be advisable to check whether the trends found in the studies by Muno and collaborators still 
hold for a larger sample. The energy dependence of the burst oscillations of the accreting millisecond pulsar J1808 also remains 
to be analysed. Our study also highlights the need to look at variations in properties over the course of the outburst; both 
accretion-powered and burst oscillation properties were found to vary fo r J1814. The technique of folding t ogether multiple 
bursts, for example (used both in this paper and in the previous study bv lMuno. Ozel fc Chakrabartv I i2003i) 'l. could obscure 
interesting variations. 

The drop in fractional amplitude with energy observed in the burst oscillations of J1814 is inconsistent with the predictions 
of both surface mode and one-temperature hot spot models. If these models are correct then there must be substantial changes 
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Figure 8. En ergy dependence of fractional ampli t ude fo r four example bursts. Burst numbers accord with those in TablelTlof this paper 
and Table 1 of lWatts. Strohmayer Markwardtl 120051) . For bursts 10 and 17 a linear relation between fractional amplitude and energy 
is a good fit, although the gradients are very different. For bursts 11 and 20 the amplitude is relatively constant between the 2-5 and 
5-10 keV bands, dropping only in the 10-20 keV band. 



in the accretion process, triggered by the burst, that cause the accretion pulsations to mask this property of the thermonuclear 
pulsations. If we assume that the accretion process does not alter significantly during the burst then we must seek an alternative 
theory. One possibility is that the hot spot possesses a temperature gradient. In Section [3.31 we examined the magnitude of 
the gradient that might be required to explain the observations. It is substantial: a rise of ~ 10% from center to edge. How 
such a gradient might be generated and whether it can meet other constraints posed by, for example, the pulse profile, are 
both interesting questions for future work. Other processes such as beaming due to a corona, however, may have to be invoked 
to explain the observations. 
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bursts. No attempt has been made to correct for the accretion contribution (see discussion in the text). 
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Figure 11. Phases calculated for the folded set of bursts, relative to the best fitting constant phase. The top panel shows a higher energy 
resolution than the lower panel. Both sets of data are consistent with there being no phase shift between the different energy bands. 
Note that 0.01 cycles corresponds to 31.8/^s. 
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Figure 12. Rotational modulation amplitude as a function of energy for two circular hot spots; one with uniform temperature, and one 
with a linear variation in temperature across the spot. In these calculations we used a spot size (half angle) of 138°, a location of the 
spot center of 80°, and an observer viewing angle of 45°. The temperature at the center of the spot is 2.9 keV, whilst the temperature of 
the rest of the star is 0.8 keV (although there is little change if this is lower). The upper trace shows the result for a uniform temperature 
hot spot. The lower trace shows the result for A = 0.1, so that the outer boundary of the spot is at 3.2 keV. 



